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We represent a new microscopic approach that allows studying the propagation properties of microwaves in a Bose- 
Einstein condensate of alkali atoms. It is assumed that the frequency of signal is tuned up to the transition between 
hyperfine ground state levels of such atoms. Pulse slowing conditions dependence on the system parameters is found. It 
is shown that the slowed signal can propagate in mentioned system with rather small energy loss. Such phenomenon is 
also studied in case of hyperfine levels Zeeman splitting. A possibility of ultra-slow microwaves observing in a condensed 
gas of cesium atoms is discussed. 

Keywords: Green functions, alkali atoms, hyperfine structure, Bose-Einstein condensate, ultra-slow electromagnetic 
waves. 

PACS: 05.30.-d; 03.75.Hh; 42.25.Bs. 



1 Introduction 

In 1999 the unique opportunity to observe the ultra- 
slow pulses of the optical region of spectra in a Bose- 
Einstein condensate (BEC) of sodium atoms was shown 
PP. In this and subsequent experiments laser pulses 
were tuned up to dipole-allowed transitions. 

Naturally, a question arises about the possibility 
of using another levels for observing electromagnetic 
pulses slowing in BEC of alkali atoms. In our opinion, 
it seems to be convenient to use microwave pulses tuned 
up to the transition between ground state hyperfine 
structure levels of hydrogen-like atoms. It is easy to 
see that such levels are rather stable. This fact, in par- 
ticular, explains the universe abundance of microwaves 
with "famous" 21 cm wavelength, which correspond to 
the transitions between hyperfine states of the atomic 
hydrogen. Moreover, there are many trapping-related 
experiments with ultra-cold atoms, which are prepared 
in different hyperfine states (see, e.g., experiments for 
the multi-component BEC observing [5]). 

However, as it is easy to see, alkali atoms in the 
ground state do not have a dipole moment. As a result, 
theoretical descriptions (see, e.g., ref. ,3j), in which 
atoms are considered as dipoles, become no longer con- 
venient. For our opinion, the most appropriate way to 
go out from this situation is using a microscopic ap- 
proach. Such an approach, for example, was developed 
in ref. [4]. It is based on the approximate formulation 
of the second quantization method and takes into ac- 
count the presence of bound states of particles. As it is 
shown below, such method is universal in some respect 
and allows to describe the response of the system not 
only to microwave range perturbation (the signal tuned 
up to the transitions between hyperfine ground state 



levels), but also to the optical range perturbation (the 
pulse tuned up to the dipole-allowed transitions) [5]. 

Let us briefly state out some basic principles of the 
suggested approach. 

2 Theoretical basis 

As one can suggest, alkali atoms are similar to the hy- 
drogen atom by the internal structure, and can be con- 
sidered as bound states of two " elementary" particles 
(valence electron and atomic core). From the stand- 
point of the quantum mechanics, it is also known that 
a bound state energy can lie not only in the discrete 
spectrum region, but also in the continuum one. This 
fact was taken into account in ref. [4j in construction of 
novel formulation of the second quantization method in 
the presence of bound states (atoms). It was also used 
in ref. [5] in construction of microscopic theory of the 
hydrogen-like low temperature plasma response to the 
external electromagnetic field. 

Thus, to construct a microscopic approach one needs 
to consider a system, that consists of two kinds of free 
fermions (ions and electrons) and their bound states 
(alkali atoms), as it was done in ref. [I]. But, to sim- 
plify the further description, we have to note the follow- 
ing. In recent theoretical investigations [6] of hydrogen- 
like low temperature plasma in an equilibrium state it 
was shown that the density of free fermions is exponen- 
tially small in comparison with the density of atoms. 

Therefore, to build microscopic approach describ- 
ing electrodynamic processes at low temperatures one 
can use the mentioned second quantization method 
with preserving the quantities that correspond only 
to bound states (atoms) contribution. According to 
ref. [4], the system Hamiltonian can be written in the 



following form: 

H(t) = H + H int + V(t), H =H ph + H p , (1) 

where Ti. p h is the Hamiltonian for free photons and 7ii n t 
is the Hamiltonian of interaction between atoms (we 
neglect Hamiltonians 7i p h and 7Yi n t below, but for the 
explicit form see also ref. [4]). 

The operator 7i p in eq. (fTJ) is the Hamiltonian for 
free particles (atoms) 
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(2) 



where (X), fj a (X.) are the creation and annihilation 
operators of hydrogen-like (alkali) atoms with the set 
of quantum numbers a at the point X; e a is the en- 
ergy of an atom in the state characterized by quantum 
numbers a; to is the atomic mass. 

And, finally, the operator V(t) in eq. |T]) represents 
the Hamiltonian that describes the interaction of par- 
ticles with the electromagnetic field 



V(t) 



dxA (e) (x,i)j(x) + / dxy>^(x,t)S-(x), 



where (x, t) and ip^ (x, t) are the vector and scalar 
potentials of the external electromagnetic field, respec- 
tively, and operators j(x) and <r(x) are the charge and 
current density operators, respectively: 
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Here V is the system volume. Note that the charge 
and current density matrix elements in eq. ([3|) can be 
expressed (see ref. [4] ) in terms of atom wave functions 
Pa(x): 
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where e is the electron charge absolute value, m p and 
to c are the atomic core and electron mass, respectively 
(TO = TO p + to c ). 



3 System response to the ex- 
ternal electromagnetic field in 
the framework of Green func- 
tions formalism 

Non-relativistic equations of quantum electrodynam- 
ics, which were found in ref. [4] on the basis of Hamil- 
tonians Ilj-Jll) and which, in turn, were used in ref. [5], 
allow to study the linear response of the system to a 
perturbation by the external electromagnetic field [7]. 
To find out such a response it is most convenient to 
use the Green functions formalism (as it was done in 
ref. [5]). In the framework of this formalism the scalar 
(charge) Green function can be defined as: 

G (+) (x,t) = -i6(t)Spw[a(x,t),fr(0)], 

where 6(t) is the Heaviside function, w is the Hibbs 
distribution operator and charge density operators (see 
definition (|3|)) must be taken in the Heisenberg repre- 
sentation. Neglecting the interaction between particles 
(as it can be done for dilute gases), one can get the 
expression for the scalar Green function Fourier trans- 
form (see ref. [5]): 
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where e Q (p) = e a +p 2 /2m 7 f a (p) is the bosonic distri- 
bution function of the ideal gas of hydrogen-like (alkali) 
atoms 

/ a (p) - {exp[(£ Q (p) - fi a )/T] - I}" 1 . 

Here [i a is the chemical potential of atoms in the state 
a, T is the temperature of the gas that is taken in 
energy units. 

Using the developed theory, it is not difficult to find 
the permittivity of such gas at low temperatures. From 
eq. ([5]) one gets (see also ref. [5]): 
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As it is known, at extremely low temperatures a 
Bose-Einstein condensate (BEC) of alkali atoms can 
be formed. At temperatures much lower the critical 
point temperature T <C To the atomic distribution 
functions f a (p) are proportional to the Dirac delta- 
function 6(p). Note also that at temperatures T < Tq 
the chemical potential [i a does not depend on tempera- 
ture and must be set equal to energy of the lowest level 
of atoms, which form BEC (see refs. [6l [8] ) . Therefore, 



after integration of eq. ^ over momentum p the ex- 
pression for the permittivity of the studied gas in BEC 
state (T — > 0) takes the form: 
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where v a is the density of condensed atoms in the quan- 
tum state a, £k = k 2 /2M, and the quantities o a p (k) 
are still defined by the formula ([4]). Note that due to 
damping processes in real systems we also introduce 
the linewidth j a /3, concerned with the transition prob- 
ability between the state a and state (3. As it easy to 
see, in eq. (O at frequencies close to the energy inter- 
val Ae a (3 (Ae a (3 = e a — ep) some peculiarities appear. 
In fact, such a behavior must have a strong impact on 
the dispersion characteristics of the gas and can be an 
underlying condition for pulses slowing. 

To study the propagation properties of the signal one 
should add the dispersion relation for free waves, that 
can spread in the studied system: 



-e(k, w)/i(k, uj) - k 2 = 0, 



(8) 



where ^(k, u) is the magnetic permeability that can be 
found analogically in the framework of Green functions 
formalism (see in that case ref. [5]). 

4 Microwaves slowing caused by 
hyperfine structure levels 

In the previous section the expression for the permit- 
tivity of the system in BEC state (see eq. (J7])) was 
represented. For the system with the frequency of the 
external field tuned up to the difference between two 
defined levels (marked below by subscripts 1 and 2) it 
can be written in a more suitable form: 
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Here gj is the degeneracy order over total spin of j 
level (j = 1, 2), Soj = ui — As2i is the laser detuning 
(\Suj\ <C A£2i), 7 = 712 is linewidth related to the tran- 
sition probability from the upper to lower state. Note 
that we also neglected the term (see eq. ([7])), the 
substantiation of such operation is discussed below. 

To study the propagation properties it is more conve- 
nient to turn to the refractive index and the damping 
factor quantities. To do it we set the magnetic per- 
meability in the dispersion relation ([5]) close to unity, 
/i(k,cj) = 1. In this case the refractive index n(k,w) 
and the damping factor x(k,w) can be expressed in 



terms of the real and imaginary part of the permit- 
tivity (e' and e", respectively), which in turn can be 
derived from eq. ©: 
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Now one can find the dependence of group velocity 
on system parameters. As it is known, the group ve- 
locity of a propagating pulse can be defined as: 



n + uj(dn/duj) 



Thus, after some mathematical transformations, in 
case of energy dissipation smallness and strong disper- 
sion, one gets: 
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This expression gives the opportunity to study the de- 
pendence of the group velocity not only on detuning Suj 
and linewidth 7, but also on characteristic properties 
of the system under consideration that have an impact 
on the parameter a value (see eq. (fTTj) '). 

Now let us obtain conditions, in which the ultra-slow 
microwaves phenomenon for two-level system in BEC 
state can be observed. To this end we proceed to the 
limit Suj — » 0. In this case, according to eq. (Tit))) , the 
real and imaginary parts of permittivity tend to the 
limits: 
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Hence, the dissipation smallness condition (e' <C |e"|) 
can be written as follows: 

M « 
— < 1. 

7 

Note that it is necessary also to add the slowing down 
(strong dispersion) condition, which in this case (see 
cq. (H2)) can be written as: 
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For the defined system with a fixed energy structure 
(Ae, 7 = const) the only parameter that can be varied 
is the occupation difference (y\ —V2), which is included 
in the parameter a (see eq. (fTTj)). Thus, basing on 
these relations we get the expression that characterizes 
the region where the mentioned phenomenon can be 
observed: 
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(13) 



By considering the example of the hyperfine levels of 
the ground state for cesium atoms let us demonstrate 
that such region can exist. The choice of such levels, 
as it is mentioned above, is stimulated by their stabil- 
ity and pumping capability. Note also that for such 
levels the dipole transitions are forbidden, thus transi- 
tions come from the higher order effects that result in 
extremely small values of linewidths. It is shown be- 
low that such fact gives the opportunity for a signal to 
propagate with a small loss of energy. 

It should be mentioned that the description can be 
extended to other hydrogen-like atoms and other type 
of levels analogically. For example, it can be used for a 
description of experiments with an ultra-slow light in 
BEC of sodium atoms pQ, in which the dipole-exited 
states were used for a pulse slowing. In this case to 
get numerical results from the developed approach one 
can restrict calculations of the charge density matrix 
element <7i 2 (see eq. 2]) to the first order approximation 
over (ky) <C 1: 

cr$(k) ~ «kdi2, 

with the dipole moment d 12 that corresponds to the 
dipole transition 1 — > 2: 

di2 = e J ydy^i(y)<y2 2 (y). 

But, it is well known that alkali metals ( 133 Cs in par- 
ticular) in the ground state do not have a dipole mo- 
ment d, thus to describe the pulse slowing on dipole- 
forbidden transitions the charge density matrix ele- 
ment a 12 (see definition (j4])) must be expanded to the 
second order over (ky) <C 1. As a result, one gets: 

«f(*n,) 2 , (14) 

where vq is the atomic radius (for cesium ground state 
7- w 2,6 x 10 -8 cm [§]), e is the electron charge. Tak- 
ing gi — 7, 32 = 9 (gj = 2Fj + 1, Fj is the total spin of 
an atom in j hyperfine state, j = 1,2), k = (Ae 21 /c), 
where Ae 2 i ~ 3, 8 x 10 -5 eV (microwaves with fre- 
quency 9,1926 GHz), the linewidth 7 w 3, 8 x 10" 21 eV, 
corresponding to the anticipated accuracy (10 -16 ) in 
"cesium fountain clock" experiments [10], and basing 
on the expressions (fTTj) . (fT3| one can find the region 
for the occupation difference, in which the pulse slow- 
ing phenomenon can be observed: 

1CT 3 cm" 3 < \v! - v 2 \ < 3 x 10 13 cm" 3 . (15) 

From this inequality one can conclude that the effect 
becomes greater with the density difference increasing 
until it reaches the upper limit of the expression (|15l) , 
when damping effects prevail in the system. We should 
stress that the region of densities (fT5|) looks convenient 
from the standpoint of experiments for BEC regime 
(see, e.g., the experiment [llj . in which cesium atoms 
with the density v — 7 x 10 10 cm -3 were used). 



We should also note the following. As it easy to see 
from eq. (fT2|) in the limit 5u — > 0, the sign of the group 
velocity v g depends directly on the sign of the quantity 
a, that in turn depends on the sign of the difference 
(y\ — 1/2). In other words, it depends on whether the 
population is normal or inverse. 

In case of normal population (i>i > v-i) the group 
velocity of signal is negative. It is traditionally consid- 
ered that the group velocity for the transparent matter 
is positive. But here one can conclude that due to the 
relation (|13|) . the signal can propagate in the system 
with rather small dissipation (i.e., in fact, the mat- 
ter is transparent) and rather slow velocity. Let us 
note that an observing of electromagnetic pulses with 
negative group velocity is not so abnormal. The exis- 
tence of such kind of phenomena for physical systems 
when the wave frequency is close to atomic (or molec- 
ular) resonances was pointed out and studied in many 
works (both theoretical, see e.g. [12], and experimen- 
tal OH!]). In case of inverse population (y\ < v%) 
more " normal" situation takes place because the group 
velocity of the slowed pulse is positive. 

We stress that such rather unusual phenomenon oc- 
curs due to the unique property of hyperfine splittcd 
ground state levels. One can show that for a two-level 
system with allowed dipole transitions the pulse will 
not propagate due to a large absorbtion. The most ob- 
vious way to go out from this situation is the second 
(coupling) laser usage that leads to the electromagneti- 
cally induced transparency (EIT, a detailed description 
see in ref. [15]). 

Now, let us say a few words about the quantity £k, 
which was neglected in deriving the equation ([9]). If 
we assume k = (Ae 2 i/c), one can find that, e.g., for 
cesium, m 3, 5 x 10~ 30 e V. So, even at the point 
Slo = it is small in comparison with the linewidth 
7 = 3, 8 x 10~ 21 eV. Thus, we proved that the used 
approximation is correct. 

5 Pulse slowing concerned with 
Zeeman splitted hyperfine 
structure levels 

In the previous section the observing possibility for 
ultra-slow microwaves caused by hyperfine structure 
levels was shown. There we assumed that such lev- 
els are degenerate, i.e. the external magnetic field is 
absent. But the presence of an external magnetic field 
(such situation occurs in most of experiments) results 
in more complicated picture. Each of hyperfine struc- 
ture levels splits to additional components with differ- 
ent total spin projection mp (e.g. for cesium atoms 
see fig. [l]). As it is well known, in a weak external 
magnetic field the energy difference between such com- 
ponents is proportional to the magnetic field intensity 
(the Zeeman effect). 
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Figure 1: Ground state energy structure of cesium atom 
in the external magnetic field. The first sign in brackets 
corresponds to the total atomic spin F, the second one to 
its projection tuf- For good visibility and compactness not 
all sublevels are showed. 



If the components of the upper hyperfine structure 
level are kept away from each other (the linewidth of 
levels is much less than the energy difference between 
levels of multiplct 5e, as shown in fig. []}, the signal 
tuned up exactly to the transition between the occu- 
pied lower state and one of the upper states (narrow 
black arrowed lines in fig. []} can be slowed down, as 
it is pointed out in the previous section. The only 
difference is that the states are no more degenerated 
(.9i = ff2 = 1) and the linewidth jj in eq. {7J corre- 
sponds to the transition probability from the chosen 
component with the projection j = (-F', ...,F') to the 
occupied lower state. Taking all other values the same 
as when deriving the inequality (fT5|) in the presence of 
the external magnetic field one can find: 



10" 
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(16) 



In the experiment [16] for the BEC regime cesium 
atoms with a peak density of 1,3 x 10 13 cm" 3 were 
kept in a trap. Thus, according to the inequality (|16[) . 
one can conclude that in such system the ultra-slow 
microwaves phenomenon could be observed. Moreover, 
direct calculations show that the group velocity of the 
signal propagating in the mentioned system is close 
to 9 x 10~ 6 m/s. It means that the pulse tuned up 
in such way covers the distance of 1 millimeter in the 
condensed cesium vapour nearly in 2 minutes! 

Let us note that to observe such a kind of phe- 
nomenon one needs to use the microwave signal with 
the linewidth 7 S much less than the level linewidth 



7s < 7i- 



(17) 



Therefore, one can meet essential experimental difficul- 
ties because, as it was mentioned earlier, the linewidth 



concerned with the transition between hyperfine struc- 
ture levels is sufficiently small. 

However, in the presence of an external magnetic 
field there are more regions where slowing down phe- 
nomenon can be observed. For example, the microwave 
signal, detuned relatively to the transitions between 
two neighbour states of the upper hyperfine multiplet 
and occupied lower state (wide grey arrow line on fig. [I} 
can also be slowed down to sufficiently small values. In 
that case the group velocity to a greater extent depends 
on the magnetic field intensity and the condition (|17|) 
is not necessary. 

On the basis of the developed approach it is not dif- 
ficult also to describe such variant of slowing, taking 
into account the Zeeman splitting. If we assume that 
all atoms are condensed in lower (occupied) state with 
the density v and energy e , we get (see eq. 0): 
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where oi2(k) can be found from the expression ()14|1 
and Aej — (£\F',m ,> ~~ £ o)- Now, using the derived 
equation, one can study the dependence of the refrac- 
tive index on the frequency of microwave signal and its 
slowing down conditions. 

Let us demonstrate it also by the cesium atoms ex- 
ample. The dependence of the refractive index in case 
of the external magnetic field presence has a rather 
complicated form due to a large number of splitted 
levels (9 levels in the upper state with F = 4), thus it 
is more convenient to show such dependence in a fig- 
ure. In fig. [2] one can see that in central part (Su) = 
that corresponds to the wide grey arrowed line in fig. [I} 
the refractive index has a slope with steepness depend- 
ing on the distance Ss between the levels |4, —4) and 
|4, — 3). It means that the microwave signal, which 
is detuned in such way, can also be slowed. The di- 
rect calculations show that, e.g. for the ideal gas of 
cesium atoms in BEC state (y = 1,3 x 10 13 cm~ 3 , 
7j = 3, 8 x 10 ~ 21 eV) with the energy of splitting 
6s = 3,8 x 1CT 18 eV (Se/jj = 10 3 , such as showed 
on fig. H]), the pulse can be slowed down to 1,07 m/s. 



6 Conclusion 

Thus, by means of the microscopic approach we studied 
the linear response of an ideal atomic gas to disturbing 
effect of an external electromagnetic field. Our ap- 
proach was based on a novel formulation of the second 
quantization method in the presence of bound states of 
particles [3] . The use of such an approach allowed us to 
obtain expressions for the dielectric permittivity of di- 
lute gases of alkali atoms in BEC state. The existence 



n(co) 



-1.5 



1.001 



1.0005 



0.9995 



0.999 



Sco/Se 



Figure 2: Refractive index dependence for cesium atoms 
in BEC state in the region of frequencies close to the tran- 
sitions between hyperfine splitted levels. Left and right 
narrow peaks correspond to the transitions |3, 3) — + |4, —4} 
and |3, 3} — > [4, —3), respectively. Note that with frequency 
increasing (to the right side from the figure edge) the be- 
havior qualitatively repeats. The used parameter for the 



Zeeman splitting energy is Se/^j 



1(T 



of resonance frequencies in expression for permittivity 
was found. 

Our approach raised a possibility to study propaga- 
tion properties of the microwave signal, tuned up to 
the transition between two hyperfine ground state lev- 
els of alkali atoms that were considered in BEC state. 
In contrast to dipole-allowed transitions, it was demon- 
strated that the pulse could propagate in such system 
with rather small energy loss. Due to this fact we intro- 
duced the group velocity concept. The slowing down 
conditions for the signal that propagated in BEC (at 
the limit of zero temperatures) were studied. More- 
over, we revealed the dependence of the group velocity 
sign on the level's population difference. Thus, we sug- 
gested that in some cases in a Bose-Einstein condensate 
could propagate the weak-damping microwaves with a 
negative group velocity. 

In case of the external static magnetic field (we took 
into account the Zeeman splitting) the ultra-slow mi- 
crowaves phenomenon was also studied. Considering 
the example of cesium atoms vapour it was shown that 
in some conditions the pulse could be slowed down to 
extremely small values. 
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